ON CENTRE SUBSPACE BEHAVIOUR 
IN THIN FILM EQUATIONS 



V A GALAKTIONOV t AND P J HARWINJ 

Abstract. The large-time behaviour of weak nonnegative solutions of the thin film 
equation (TFE) with absorption 

m = -v ■ (M n vAu) - 

with parameters n S (0,3) and p > 1, is studied. The standard free-boundary prob- 
lem with zero-height, zero contact angle, and zero-flux conditions at the interface and 
bounded compactly supported initial data is considered. It is shown that there exists 
the critical absorption exponent 

Po = 1 + n + jj 

such that, for p = pq, the asymptotic behaviour of solutions u(x,t) for t 3> 1 is repre- 
sented by the well-known source-type solution of the pure TFE absorption, 

u s (x,t) = t~ /3N F(y), y = x/t , with the exponent (3 = j^jj, 

which is perturbed by a couple of lnt-factors. For n = 1, this behaviour is associated 
with the centre subspace for the rescaled linearized thin film operator and is given by 

u(x,t) ~ (tint)-"* FCs/*" (In with = 

where F(y) — 8nv+2T71v+4l ( n * — l^l 2 ) 2 an< ^ ^ ne constant a* > depends on dimension 
TV only. The 2mth-order generalization of such TFEs with critical absorption is consid- 
ered and some local and asymptotic features of changing sign similarity solutions of the 
Cauchy problem are described. 

Our study is motivated by the phenomenon of logarithmically perturbed source-type 
behaviour for the second-order porous medium equation with critical absorption 

u t = V • (u n Vu) - up inR N xR+, p = 1 + n + n > 0, 

which has been known since the 1980s. 



I. INTRODUCTION: THE MODEL, MOTIVATION, AND RESULTS 

Our goal is to describe some unusual asymptotic phenomena for higher-order quasilin- 
ear degenerate parabolic equations, in which the nonlinear interaction between operators 
involved deforms the scaling-invariant structure of solutions for large times. These del- 
icate cases of asymptotic phenomena, such as logarithmic perturbations of fundamental 
or source-type solutions, have been known since the 1980s for quasilinear second-order 
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reaction-diffusion equations. For semilinear higher-order parabolic equations, those phe- 
nomena can be detected by using spectral theory of non self-adjoint operators and semi- 
group approaches. For quasilinear models, similar asymptotic patterns were unknown. 

In the present paper, we introduce a new quasilinear parabolic model by adding to the 
standard thin film operator an extra absorption term. This creates a non- conservative 
evolution PDE, which enjoys a variety of logarithmically perturbed non-scaling asymp- 
totics in both free-boundary and the Cauchy problem. We then fix several similarities 
with simpler second-order diffusion-absorption models. 

We begin with some physical motivation of such models. 

1.1. On general thin film models: a class of conservative and non-conservative 

PDEs. For a long time, modern thin film theory and application dealt with rather com- 
plicated nonlinear models. Typically, such models include the principal quasilinear fourth- 
order operator and several lower-order terms. For instance, the Benney equation (1966) 
describes the nonlinear dynamics of the interface of 2D liquid films flowing on a fixed 
inclined plane [2J, 

(1.1) u t + 2 -f(u 3 ) x + e[(^ffV - cot 6 u 3 )u x + Zu 3 u xxx ] x = 0, 

where Re is the unit-order Reynolds number of the flow driven by gravity, a is the rescaled 
Weber number (related to surface tension a), 9 is the angle of plane inclination to the 
horizontal, and £ = | < 1, with d being the average thickness of the film and A the 
wavelength of the characteristic interfacial disturbances. See [39J. 

TFEs can include non-power nonlinearities. For instance, in the multi-dimensional 
geometry, a typical example is 

(1.2) Ut + V- [(-Gu 3 + 2P B { ^ uY )Vu] +SV- (u*VAu) = 

that describes, in the dimensionless form, the dynamics of a film in R 3 subject to the 
actions of thermocapillary, capillary, and gravity forces. Here, G, M, P, B, and S are 
the gravity, Marangoni, Prandtl, Biot, and inverse capillary numbers respectively. On 
Marangoni instability in such TFE models, see [37] . 

The above conservative PDEs preserve the finite mass of thin films. Non-conservative 
TFEs occur for evaporating/condensing films and via other effects, J38j[27]. Actually, the 
first study of the vapor thrust effects in the Rayleigh-Taylor instability of an evaporating 
liquid- vapor interface above a hot horizontal wall was performed by Bankoff in 1961. 
His stability analysis in 1971 of an evaporating thin liquid film on a hot inclined wall 
extended earlier results of Yih (1955, 1963) and Benjamin (1957). The history and detailed 
derivation of models of (a) evaporating thin film and (b) a condensing thin film, can be 
found in [381 PP- 946-949]. A typical TFE of that type in ID is as follows [381 P- 949]: 

Here, the six terms represent, respectively, the rate of volumetric accumulation, the mass 
loss, the stabilization capillary, van der Waals, vapor thrust, and thermocapillary effects. 
In the second absorption-like term, E is the scaled evaporation number and K is the 
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scaled intefacial thermal resistance that physically represents a temperature jump from 
the liquid surface temperature to the uniform temperature of the saturated vapor. D is 
a unit-order scaled ratio between the vapor and liquid densities. 

Another origin of non- conservative TFEs with more complicated non-divergent oper- 
ators is the study of flows on a rotating disc (centrifugal spinning as an efficient mean 
of coating planar solids with thin films). This gives extra absorption-like, spatially non- 
autonomous terms in the equations written in radial geometry, e.g., [381 P- 955] 

+ [rV + eRe^ErW - ^rV)] 

+ §{Re(fr 3 « 6 - r^u 3 )u r + r±u 3 [\.{ru r ) r ] r ) r = 0. 

Here E is again the evaporation number, F is the Froude number, and e = ^ is a small 
parameter. Observe a rather complicated combination of various absorption and reaction- 
like non-divergent terms (with different nonlinear powers u 3 , u A , and u 7 ) in the first line 
of equation (11.41) . Various exact solutions of non-conservative TFEs can be found in [2U 
Ch. 3], where more references and a survey on TFE theory are given. 

Modern nonlinear parabolic theory and application to thin film models demand better 
understanding of interaction of various nonlinear terms and operators of different orders 
that can create rather complicated spatio-temporal patterns and dissipative structures. 
We chose one particular but special case of centre subspace behaviour that will be shown 
to have rather robust mathematical significance. 

1.2. Basic limit model: the TFE with absorption. We study the large-time asymp- 
totic behaviour of nonnegative solutions of the thin film equation (TFE) with absorption 
(for convenience, it is written for solutions of changing sign to be studied also) 

(1.5) u t = -V- (\u\ n VAu) - [u^u, 

where n > and p > 1 are fixed exponents. Here we use the simplest second term which 
is not a differential operator but is represented by just a power function. Our main goal 
is to justify that in the critical case 

(1-6) p = l + n + ± 

various solutions of (11. 5p exhibit a complicated asymptotic behaviour with some logarith- 
mic corrections hit for t > 1. 

We have chosen the non- conservative equation (I1.5P for simplicity and for better pre- 
sentation of our mathematical tools. We claim that similar phenomena are quite general 
and appear also in various conservative models. Actually, the logarithmic correction 
~ (lni) -1 / 7 in the behaviour for large enough t was rigorously observed [2E] for the 
relaxed conservative thin film model consisting of two operators, 

(1.7) u t + {u 3 u xxx ) x + {u n u xxx ) x = 0, with 0<n<3 (u > 0), 

where the first term with u 3 corresponds to Reynolds'equation from lubrication theory. It 
was shown that, for concentrated enough initial data, in a certain intermediate time-range, 
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the propagation rate is as follows: 

(1.8) meas{u(x,t) > 1} ~ 

where the usual scaling-invariant factor t$ is associated with a standard dimensional anal- 
ysis. Here, the log-correction is a result of a delicate interaction of two scaling invariant 
operators in (jl.7p . We believe that (11.81) . proved in rigorously, can be put into a 
framework of a centre manifold calculus (though a justification can be extremely hard). 
Log-corrections were observed for the limit stable Cahn-Hilliard equation [T91 Sect. 5.4] 

(1.9) u t = -A 2 u + A(\u\ p ~ l u), with p=l + ^. 

For the semilinear case n = in the TFE (11.51) . such logarithmically perturbed asymptotic 
are also well known and admit a rigorous mathematical treatment, |21j . 

Thus, we consider for ( II. 5p the standard free-boundary problem (FBP) with zero-height, 
zero contact angle, and zero-flux (conservation of mass) conditions 

(1.10) u = Vu = v ■ (u n VAu) = 

at the singularity surface (interface) To[m], which is the lateral boundary of suppw with 
the outward unit normal v. Bounded, smooth and compactly supported initial data 

(1.11) u[x, 0) = u Q (x) in r [u] n {t = 0} 

are added to complete a suitable functional setting of the FBP. As usual, we assume that 
these three free-boundary conditions give a correctly specified problem for the fourth-order 
parabolic equation, at least for sufficiently smooth and bell-shaped initial data, e.g., in 
the radial setting. 

Returning to basics of thin film theory, earlier references on derivation of the pure 
fourth-order TFE 

(1.12) u t = -V- (\u\ n VAu) 

and related models can be found in [2"8l l4Tj . where first analysis of some self-similar 
solutions for n = 1 was performed. Source-type similarity solutions of (11.121) for arbitrary 
n were studied in [7] for N — 1 and in [20] for the equation in M . More information on 
similarity and other solutions can be found in [5^14"! Illj. In general, the TFEs are known to 
admit non-negative solutions constructed by special "singular" parabolic approximations 
of the degenerate nonlinear coefficients; see the pioneering paper [3J, various extensions 
in [291 EH EE E21 Si] an d the references therein. In what follows we study the asymptotic 
behaviour of sufficiently "strong" weak solutions of the TFEs, which satisfy necessary 
regularity and other assumptions; see also the survey paper [TJ. Notice that regularity 
theory for the TFEs is not fully developed, especially in the non-radial iV-dimensional 
geometry and for solutions of changing sign, so we will need to impose extra formal 
requirements, which are necessary for justifying our asymptotic approaches. 

Let us mention other well-established and related conservative thin film models with 
extra lower-order terms describing the dynamics of thin films of viscous fluids in the 
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presence of two competing forces; see [9]. For N — 1, typical quasilinear TFEs are 

(1.13) u t = -{uu xxx + u 3 u x ) x (u>0), 
and the general equation with power nonhnearities is 

(1.14) u t = -{u n u xxx ) x -{u m u x ) x (u>0). 

We refer to papers [T71 HE] and the book [2U Ch. 3] as sources of a large amount of further 
references and results of TFE theory and application. 

In addition, our extra motivation of the TFEs model like (jl.5p is mathematical and is 
associated with the previous investigations of the quasilinear diffusion-absorption PDEs. 

1.3. A mathematical motivation: the PME with critical absorption. Second- 
order quasilinear parabolic equations with absorption are well known in combustion the- 
ory. A key model is the porous medium equation (PME) with absorption 

(1.15) u t = V • (u n Vu) - u p in xl + (u > 0), 

where n > and p are fixed exponents. A special interest to such equations was mo- 
tivated by localised similarity solutions introduced by L.K. Martinson and K.B. Pavlov 
at the beginning of the 1970s. Mathematical theory of such PDEs was developed by 
A.S. Kalashnikov a few years later; see his survey [30] for the full history. Besides new 
phenomena of localization and interface propagation, for more than twenty years, the 
PME with absorption (11.151) became a crucial model for determining various asymptotic 
patterns, which can occur for large times t ^> 1 or close to finite-time extinction as t —>■ T~ 
(for p < 1). For (j!.15p . there are a few parameter ranges with different asymptotics, 

p>p = l + n + jf, p = p , l + n<p<p , p=l + n, 
1 < p < 1 + n, p — 1, 1 — n < p < 1, p—l — n, p < 1 — n, 

etc.; see references and details in [251 Ch. 5,6]. 

The most interesting and unusual transitional behaviour for f II . 1 5 j) occurs at the first 
critical (or Fujita) absorption exponent 

(1.16) p = l + n + l. 

In this case (see details and references in [25j p. 83]), the asymptotic behaviour as t — > oo 
of nonnegative compactly supported solutions of (11.151) is described by the logarithmically 
perturbed source-type solution of the pure PME, 

(1.17) u(x,t) = (tlnt)-^[F(x/^(lnt)-^/ 2 ) + o(l)], where /3 = 

Without the logarithmic factors and the o(l)-term, the right-hand side is indeed the 
famous Zel'dovich-Kompaneetz-Barenblatt (ZKB) similarity source-type solution of the 
pure PME Ut — V • (ii n Vii), which has the form 

(1.18) u s (x,t) =t^ N F(y), y = x/iP i with F(y) = [f (a 2 - \y\ 2 ) + ] - , 
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where a > is an arbitrary scaling parameter. This explicit solution dates back to the 
1950. In the class of solutions of changing sign, (I1.15P admits a countable sequence of 
critical exponents, where the patterns contain similar logarithmic time-factors, 



1.4. Outline of the paper: logarithmically perturbed patterns for the TFE with 
absorption. In Sections [2] we show that similar logarithmically perturbed source-type 
patterns exist for the TFE with absorption (11.51) . with the critical exponent (11.61) . In this 
case, the source-type solutions of the TFE (I1.12p take the form 

(1.19) u 8 (x,t)=t-P N F(y), y = x/lP, with/3 = iT ^, 



where F(y) > is a radially symmetric compactly supported solution of the PDE [20] 

(1.20) A(F) = -V ■ (F n VAF) + f3VF ■ y + (3NF = 0. 

In the case n — 1, the similarity profile for the FBP is given explicitly 

(1-21) F(y) = c (a 2 - \y\ 2 ) 2 , c = 8(JV+2 j (JV+4) , a > 0, 

and was first constructed in [41]. Figure [T] shows profiles F(y) for N = 1 in four cases 
n — j, |, | and 1. The profiles are normalised by their values at y = 0, so F(0) = 1. 

First, for n — 1, relying on the explicit representation (ll.2ip and good spectral prop- 
erties of the corresponding self-adjoint linearised rescaled operator, we show that, for 
V — Po — 2 + jr, the TFE with absorption (I1.5P admits asymptotic patterns of the 
following form: 

(1.22) u (x,t) ~ (tlnty^F^x/t^lnt)-^) (0 = -L_) . 

Here F* is a fixed rescaled profile from the family (ll.2ip with a uniquely chosen parameter 
a = a* > that depends on iV only. We also present evidence that similar logarithmic 



factors can occur for arbitrary n > but this does not lead to self-adjoint linearised 
operators and explicit mathematics. On the other hand, for the semilinear case n — 0, 
i.e., for the fourth-order parabolic equation written for solutions of changing sign 

(1.23) u t = -A 2 u-\u\ p - 1 u, 

the critical behaviour like ( 11. 22ft is known to occur at the critical exponent p = 1 + -4 
[21], which is precisely (11.61) with n = 0. In this case, the centre manifold analysis also 
uses spectral properties of a non self-adjoint linear operator studied in [131 Sect. 2]. 

In Section [3] we briefly describe the essence of the easier supercritical case p > p$. Very 
singular similarity solutions (VSSs) in the subcritical one p G (n + l,po) w iU be studied 
in a forthcoming paper. 

In Section HI we explain how the critical asymptotic behaviour occurs for the 2mth-order 
TFE with absorption 

(1.24) m = (-l) m+1 V ■ (« n VA m - x «) - u p , m> 2, 
where the critical absorption exponent is 

(1.25) p = l + n + ^, 

and again n — 1 leads to a simpler self-adjoint case. 

In Section [5] we discuss similar local and global asymptotics for the Cauchy problem 
admitting maximal regularity solutions of changing sign. 

Finally, let us note that (11.61) (and (I1.25P for equation (jl.24p ) is the critical Fujita 
exponent of the TFE with source 

Ut = -V • (|w| n VAw) + \u\ p {n>0,p>l), 

i.e., forp G (n+l,p ], all solutions with arbitrarily small initial dataw (x), where J u > 0, 
blow-up in finite time, [23] . 

2. Rescaled equation and centre subspace behaviour 

2.1. To the style of the analysis. For convenience of the Reader, we must emphasize 
from the beginning that all our final conclusions on centre subspace behaviours detected 
below are mathematically formal when we deal with the quasilinear case n > 0. The 
semilinear case n = is easier and admits a rigorous treatment by invariant manifold 
theory, |21j . It is then worth mentioning that there is no hope that such asymptotics 
can admit a reasonably simple rigorous treatment. We recall that even for the second- 
order model (11.151) with n > 0, there is no a full centre manifold justification of the main 
results that were proved by essential use of the Maximum Principle and comparison-barrier 
techniques; see [23 Ch. 4]. Some of asymptotic patterns for (11.151) of centre subspace type 
turned out to be very complicated, [22] . As we will show, the main difficulty is not a proper 
spectral theory of linearized operators (this is justified in many cases) but a justification 
of the centre subspace behaviour associated with such singular operators. On the other 
hand, we always clearly indicate the rigorous steps and split the whole approaches into a 
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sequence of standard steps. We would be very pleased if some of our formal results and 
discussions would attract attention of experts in these areas of differential equations. 
Thus, in what follows, we use by implication the following rule: 

(i) all conclusions concerning spectral and other properties of self-adjoint singular el- 
liptic and ordinary differential operators are rigorous (or can be made rigorous after 
sometimes technical manipulations; for non-self-adjoint cases we are not that certain and 
extra analysis is necessary); and 

(ii) further extensions via above spectral properties to describe the behaviour for TFEs 
close to center subspaces and various matching procedures are mathematically formal. 

2.2. Rescaled equation. We begin with rescaling the PDE (11.51) with the critical expo- 
nent (jl.6p according to the time-factors of the source-type solution (11.191) . i.e., by setting 

(2.1) u (x, t) = (l + t)^ N v(y, r), y = x/(l + tf , r = ln(l + t), 
that leads to the following autonomous rescaled equation in R N x K + : 

(2.2) v T = A(v) - v p , 

where A is the operator specified in (11.201) . We first need to check that a simple stabi- 
lization as r — > +oo to a nontrivial stationary solution in (12.21) is not possible. 

Proposition 2.1. The stationary equation 

(2.3) A(g) -g p = 

does not have a nontrivial compactly supported nonnegative solution of the FBP. 

Proof Indeed, integrating ( 12. 31) over supp g yields J g p (y) dy — 0. □ 

This means that the only bounded nonnegative equilibrium for the dynamical system 
is trivial, 

(2.4) g(y) = in R N . 

In order to detect the actual non-stationary asymptotic behaviour, we next perform 
second rescaling by introducing the as yet unknown positive function b(r), 

(2.5) v(y,r) = b(r)w(C,r), C = y/b$(r), 
to get the following perturbed equation: 

(2.6) w T = A(w) + ^Cw- bP^vjP, where Cw = \ Vw • C - w. 

2.3. Linearisation. Roughly speaking, in order to detect the asymptotic behaviour ac- 
cording to (I2.5p . we can use the estimate 

(2.7) b(r) w ||f(-,r)|| 0O —> as r — ► oo, 

so that \\w(-, t) ||oo ~ 1 for r ^> 1. On the other hand, in the radial setting, it is convenient 
to use b{r) for the scaling of the support of the solution w((, r) to get that it approaches 
the unit ball B\ as r — ► oo; see below. 
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We next perform the linearisation by setting 



(2.8) w((,t) = F(()+Y((,t), 

where F is a rescaled similarity profile from the family (ll.2ip . Then Y solves the following 
rescaled equation: 

(2.9) Y T = A'(F)Y + b T CF - b^F? + D(F) — 6 P_1 [(F + Yf - F v ] , 
where A'(F)Y is the formal Frechet derivative of A at F, 

(2.10) A'(F)Y = -V • (F n VAY) - V • {nF n ~ l YVAF) + (3VY ■ y + (3 NY, 

and D(y) is a higher-order perturbation, which is quadratic in Y — > on smooth func- 
tions. Using the elliptic equation (jl.20p for F, on integration, 

(2.11) F n VAF = (3Fy A'(F)Y = — V • (F n VAF) + (1 — n)/?V • (Y(). 

2.4. The self-adjoint case n — 1. It follows from (12.111) that n = 1 is a special case, 
where the last term vanishes. We fix a = 1 in (11.211) . so that the linearised operator is 

(2.12) A'(F)Y = -V • (FVAF) = -c V • ((1 - |C| 2 )VAF), yeB 1 = {\(\ < 1}. 
One can see that it can be written in the form 

(2.13) A'(F)Y = -c ^ [A(a(|C|)AF) + 2NAY], where o(|C|) = 1 - |C| 2 = 

so, in the topology of L 2 p (E>i), operator ( 12.121) is symmetric in C^(Bi) with good coef- 
ficients, and hence admits self-adjoint extensions. Next, using classical theory [10], we 
specify properties of its unique Friedrichs self-adjoint extension. Its domain is constructed 
by completing C^°(Bi) in the norm induced by its positive quadratic form (corresponding 
to the operator -A'(F) - c A > 0) 

(Y, W}* = c J Bi [aAYAW - (2N - 1)VY ■ VW]. 

The intersection of this Hilbert space with the domain of the maximal adjoint operator 
D((A'(F))*) = {v G L 2 p : A'(F)v G L 2 p } defines the domain of the self-adjoint extension, 
which we denote by D(A'(F)) = H p0 . In particular, for any v G H p0 , there holds 

v = on dBi, and f B a(Av) 2 < oo, 

so that H p0 C H 2 . Consider the corresponding eigenvalue problem written in the form 

(2.14) - co[A(a(|C|)AV) + 2iVA^] = in H^ . 

Since the embeddings of the corresponding functional spaces H 2 and Hq into L 2 are 
compact, [351 P- 63], we have that the spectrum a(A'(F)) is real and discrete. 

For our purposes, it suffices to detect the eigenvalues and eigenfunctions in the radial 
(ODE) setting with the single spatial variable r = \(\ > 0. The extension to the elliptic 
setting is performed by using the polar coordinates ( = (r, a) in B\, 

(2.15) A = A r + 4 r A (T , 
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where is the Laplace-Beltrami operator on the unit sphere S N_1 = dB\ in M. N . A a 
is a regular operator with a discrete spectrum in L 2 (S N ~ 1 ) (each eigenvalue repeated as 
many times as its multiplicity), 

(2.16) °{~K) = W = k{k + N - 2), k > 0}, 

and an orthonormal, complete, closed subset {Vfc(cr)} of eigenf unctions, which are homo- 
geneous harmonic k-th order polynomials restricted to S N ~ 1 . We plug (I2.15P into (I2.13p . 
where all the coefficients are radial functions, and use the separation of variables 

(2.17) 4(r,<r)=R(r)V k (<r) 

for solving the eigenvalue problem ( 12.141) . For each fixed z/^, we then arrive at a radial 
eigenvalue problem for R, which is similar to that discussed below. 

Thus we take k = in (I2.17P and consider the radially symmetric eigenvalue problem 
(12.141) . For N = 1, this problem was studied in [8], where further references are given. It 
is not difficult to check that the radial operator A'(F) has the discrete spectrum 

(2.18) a(A'(F)) = {\ k = c k(k + 2)(k + N)(k + N + 2), k = 0, 2, 4, ...}, 
where each eigenfunction ip is a (fc+2)th-order polynomial, 

(2.19) Mr) = h(r k+2 + ... + d k ) (^(1) = 0), 

where {b k } are normalization constants, so that the eigenfunction subset {ip k } is orthonor- 
mal in L 2 . In particular, 

(2.20) Mr) = b (r 2 - 1) > 0, b = -^^ (A = 0), 

where ujn = j^mwj is the volume of the unit ball in R JV . Such polynomials are complete 
and closed in typical weighted L p -spaces (a standard functional analysis result; see [13j 
Sect. 2.3] for details), and this justifies the equality in (I2.18p . Moreover, we then can 
use the eigenfunction expansion with the orthonormal eigenfunctions subset {ipk} to deal 
with solutions of the corresponding PDE. 

We next consider the rescaled equation (12.91) . which for n = 1 takes the form 

(2.21) Y T = A'(F)Y + f CF - If^F? - V • (WAY) - bv- l [(F + Y)p - F^}. 

We deal with strong radially symmetric solutions of (12.211) . where we now choose the 
normalization function h[r) in (12.51) such that 

(2.22) suppw(-,r) = B x for r > 1. 

According to equation (I2.2ip . we then need to assume that 6(r) is smooth, at least, 
for large r, though this requirement can be weaken by using a weak (integral) form of 
the PDE. We now use the converging (in l? p and in the corresponding Sobolev class) 
eigenfunction expansion of the radial solution 

(2-23) Y((,T) = Zk> Mr)MC) 

to study the corresponding centre subspace behaviour for the nonlinear operator A. This 
part of our asymptotic analysis is formal. 
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Thus substituting (12.231) into (I2.2ip and projecting onto -0o i n L 2 , we have that the first 
coefficient satisfies the following perturbed "ODE" : 

(2.24) a' = - 7l f- - l2 bf- 1 + ... , where 7l = -(CF, tfj ) p > 0, 72 = (F p , </> > P > 0. 

We omit in (12.241) the higher-order terms assuming that, for this type of behaviour, the 
non-autonomous perturbations are the leading ones. The signs of the coefficients 7^2 in 
(12.241) are essential and are easily checked by integration. 
It follows from (J22D and (12T2^]) that b(r) -> as r -> 00, so 

is not integrable at r = 00. 

Therefore, in order to have a uniformly bounded expansion coefficient ao(r), we need to 
suppose that two terms on the right-hand side of ( 12.241) annul each other asymptotically, 
so that, up to an integrable perturbation, 

(2.25) £ = -izb p - 1 + ... forr>l. 

This gives the following necessary condition for existing of such a behaviour: 

(2.26) 6(r) = 7*r^ + ... ; w here 7, = [^j^]' 1 ^ 1 ■ 

Returning to the original variables {x, t, u}, from (I2.26P we obtain the asymptotic pattern 
(11.221) . The rescaled profile F* is uniquely determined from (11.211) with a* = y^^. 

2.5. Arbitrary n G (0, |). This non self-adjoint case is more difficult. Consider the 
linearised operator (12.111) for n 7^ 1, where F > is the radial solution of the ODE (11.201) 
in Bi, see [20] for existence, uniqueness, and asymptotics. Then, for n < | [20J, 

(2.27) ^(C)~(1-KI) 2 as ICI-1-. 

Notice that there exists a one-parameter family of the solutions given by 

(2.28) F a (C) = a^F(Z), a > 0. 

Firstly, we claim that, for n 7^ 1, operator (12. lip is not symmetric in L 2 p for any positive 
weight p in B\\ see Appendix A. Secondly, we have that 

(2-29) Vo(C) = £F.(C)U = ~F - VF ■ C 

is a positive eigenfunction of (12. lip corresponding to A = 0. Observe that, with respect to 
the regularity, this eigenfunction well corresponds to that for n = 1; cf. ( 12.201) . Moreover, 
it follows that, close to the singular point |£| = 1, the radial part of (12.111) is governed by 
the singular (at dB\) higher-order operator 

(2.30) L i Y = -(8 2n Y m y, s = l-|C|, 

which is symmetric in a weighted H~ x topology (but we need a result in L 2 ). Solving 
the problem L4Y = g with natural conditions at the point s — 1, which is assumed to be 
regular, we obtain, up to compact perturbations, that 

(2.31) L 2 Y = -Y" ~ J s s~ 2n j s g = Ug =^ Y ~ L 2 1 L*g 1 
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where L 2 1 is a compact operator in L 2 . It is easy to check that the integral operator 
is bounded in L 2 for 

(2.32) n < f , 

and then Ll^L* is compact in L 2 as the product of a compact and a bounded operator. 
Therefore A'(F) has discrete spectrum in the parameter range ( 12.321) . This is not an 
optimal result since, as we have seen, the discreteness of the spectrum remains valid for 
n — 1. We use this analysis as a simple illustration of the fact that the spectrum is usually 
discrete in the non-symmetric case. 

Thus G cx(A'(F)) is an isolated eigenvalue. There is a numerical evidence that the 
spectrum is discrete for all n G (0, |); see [8], where, moreover, first six eigenvalues turned 
out to be real for N — 1. Possibly this might mean that in a special topology of sequences 
as I 2 (not related to any of L 2 ) the linearised operator can be treated as symmetric and 
self-adjoint; cf. an example in [13]. For n = in any dimension N > 1, the whole 
spectrum is proved to be real. We refer to [I3J Sect. 2], where this and other 2mth-order 
operators were studied in L 2 {R N ), i.e., for the Cauchy (not a free-boundary) problem. 

The rest of our study is formal. Once in the radial setting there exists the centre 
subspace of A'(F), we are looking for a (formal) centre subspace patterns for (12.91) 

(2.33) Y((,t) = Oo(t)MO + -- 

We assume the centre subspace dominance in the behaviour, so, as usual, other terms 
in this expansion are assumed to be negligible for r 1. Substituting (12.331) into (12. 9p . 
we next find the projection onto the corresponding adjoint eigenfunction ip^. In general, 
such an analysis becomes rigorous if we establish existence of complete, closed and bi- 
orthonormal eigenfunction subsets {ipk} and {ipl}- This is an open problem except the 
case n = 1 above and n = studied in [13]. We do not deal with the adjoint operator 
A'*(F) in this formal asymptotic analysis. The projection onto ipo yields the perturbed 
ODE (I2.24p . where the same coefficients 71,2 are determined via the standard dual L 2 
product, where ^0 is replaced by ipo- This formally leads to the same asymptotics (12. 26[) . 

The range n G [|,3). The centre subspace analysis applies also for larger n's. The 
asymptotics of similarity profiles change at n = |, where, instead of (I2.27p . 

(2.34) ^(C)~(l-|Cl) 2 [I/5|ln(l-|C|)|] § as|C|-l; 
see [20]. On the other hand, for n G (|, 3), 

(2.35) F(C)~(1-|C|)» as |C| - 1. 

This regularity is sufficient for determining the corresponding eigenfunction and the log- 
arithmic behaviour. 

For n > 3, the zero contact angle FBP does not provide us with a proper interesting 
evolution; see [20] . 
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3. On the supercritical parameter range p> p 

3.1. Exponentially perturbed dynamical system for p > p$. Let us explain what 
we expect for p > p in (11.51) . In terms of the rescaled function 



(3.1) u(x,t) = (l + t)~T+z*v(y,r), r = ln(l + t), 
the equation takes the form 

(3.2) v T = -V • (v n VAv) + ^yVv + I £_v- e-TTyr, 

where 7 = ^pr-jr > if p > Pq. Therefore the absorption term — u p in (II .5p generates 



an exponentially small perturbation in the rescaled equation (13.21) . Hence one can expect 
the convergence as r — > 00 to the rescaled similarity profile F in (I1.19P of the limit mass, 
though the passage to the limit in (I3.2p generates a number of technical difficulties. Here 
(13. 2p is formally an exponentially small perturbation of the autonomous rescaled TFE 

(3-3) v T = A(v) = -V • (v-VAv) + ^ y .Vv + ^_v. 

As usual, we gain an extra advantage in the case n — 1. 

3.2. The gradient case n — 1. It is known that, for n = N = 1, the rescaled TFE 
( 13. 3p is a gradient system, [12J. Let us construct an "approximate" Lyapunov function 
for strong solutions of the FBP in M N . Namely, we write down ( 13.21) in the form 

(3.4) v T = V • [vV(-Av + \y\ 2 )} + e"^ 

and multiply in L 2 (R N ) by (— A v )~ 1 v T , where, by definition, 

(~A v y 1 w = g if A v g = V • (vVg) = -w, 
and g — at the free boundary of v. Then integrating by parts yields the identity 

(3.5) f *;|V(-A„)-V| 2 = £[-U |V^| 2 - 5^ Jv\y\ 2 } + J, 
where J corresponds to the exponentially small term, 

(3.6) J = e" 7T / ^(-A„)- V 
Integrating (13. 5p over (0, T) yields 

J T / v\ V(-A,)- V| 2 + \ J \Vv(T)\ 2 + w ^f v(T)\y\ 2 ] <C + fi j, 

so that, if the exponential term (13. 6p J G L X (]R + ), this yields extra uniform estimates, 

v / wV(-A^)-V G L 2 (R x E + ) and Vt», y/v\y\ G L°°(R + ;L 2 ). 

Note that, obviously, ( 13.51) does not imply existence of a Lyapunov function (the non- 
autonomous PDE ( 13.41) is not a gradient system). Anyway, since ( 13.51) gives a rather 
strong estimate of v T for r 3> 1, this makes it possible to pass to the limit r — >• 00 
and establish stabilization to an equilibrium point (see the technique in [25| p. 116-117]), 
which is unique by the obvious mass-monotonicity with time of the solution. 

The symmetry of the Frechet derivative (12.121) at F looks like a certain "remnant" of 
the fact that the original PDE is a gradient system. 
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4. Centre subspace patterns for the 27?tth-order TFE 



We consider the 2mth-order TFE with absorption (11.241) with the critical absorption 
(Fujita) exponent (jl.25p . The proper setting of a standard "zero contact angle" FBP for 
the TFE includes m+1 free boundary conditions at the free boundary T = dQ(t) x M + 
is the support of u(-,t) at time t > 0), 

(4.1) u = Vu = ... = = v- VKA™-^) = 0, 

where v is the unit outward normal to dfl(t) that is assumed to be sufficiently smooth. 



4.1. Similarity solutions. The similarity solutions of the pure TFE 
(4.2) ut = (-l) m+1 V • (u n VA m - 1 «) 

take the standard form f II . 1 9 j) with 



(4.3) P 



i 



2m+nN ' 



One can see that the critical exponent (11. 25ft is precisely the one, for which the PDE 
(11.241) possesses the same group of scaling transformation. Then the rescaled profile F 
satisfies the radial restriction of the 2mth-order elliptic equation 

(4.4) A(F) = (-l) m+1 V • (i^VA^F) + (3VF ■ y + (3NF = 0. 

It seems that, for any m > 3, the questions of existence and uniqueness of a solution 
F(y) > in B\ remain open. It is clear that, for large m, a standard approach to existence 
based on a multi-parametric shooting leads to a complicated geometric analysis (though 
some general conclusions in this geometry are likely). We expect that the approach based 
on the n-branching (or a continuous homotopy connection with n — 0) via the classical 
theory [42] makes it possible to explain properties solutions, at least, for small n > 
by branching from the linear case n = (but, surely, a standard approach to smooth 
branching does not apply). For the Cauchy problem, the spectral and other properties 
of the corresponding linear operator (14.41) for n = are given in [T3], and can be used 
to clarify the behaviour for small n > 0. For the FBP (14. ip . an extra analysis of the 
linearised elliptic PDE is necessary. 

As usual, the case n = 1 provides us with the explicit solution. Writing the ODE (14.41) 
in the radial divergent form (here y is actually \y\) 

{y N ~ 1 F(A m ~ 1 F)')' = {-l) m (3{y N F)\ 

on integration we obtain A m ~ 1 F = (— l) m | (3y 2 . Integrating this linear ODE 2m-2 times 
yields the positive solution in E>\ 

(4.5) F(y) = c (l - \y\ 2 ) m , where c - 1 m 



2 (2m)U(2m+N)U ' 
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4.2. Linearised operator. We next follow the same scheme of the asymptotic analysis 
as in Section [21 Similar to (I2.10p . we introduce the linearised operator 

(4.6) A'(F)Y = (-l) m+1 V • (F^A^Y) 

+ (-l) m+1 V • (nF^YVA^F) + (3VY ■ y + (3 NY. 
Using the ODE ffl~4"j) , we have that 

■ (nF n_1 VA m_1 F) = -(3nN, (-l) m+1 nF"- 1 VA m - 1 F = -flny, 
so (14. 6p can be written in the form 

(4.7) A'(F)Y = (-l) m+1 V • (i^VA^Y) + (3N(1 - n)y ■ VY + /3N(1 - n)Y, 
and we again observe that n = 1 is a special case. 

4.3. The self-adjoint case n — 1. Plugging the profile (14.51) into (14. 7p yields the follow- 
ing symmetric form of the operator: 

(4.8) A'{F)Y = c (-l) m+1 V • ((1 - \y\ 2 ) m V A^Y) 

= c (-i) m+1 [£> m ((l - \y\ 2 )D m Y) + to(to - l)NA m - l Y], 

where D m denotes A m l 2 for to even and VA^ m_1 ^ 2 for m odd. For instance, for N — 1 
and m = 3, we have 

A'(F)Y = c (l - y 2 ) 2 [((l - y 2 )Y"')'" + QY^}. 

Having the symmetric operator (14. 8 p in C^°, we next determine its self-adjoint extensions, 
[ID] . In particular there exists the extension with discrete spectrum and polynomial 
eigenf unctions in the radial setting (the non-radial case is covered by using the spherical 
polynomials as in ( 12.17P ). The eigenvalues for the polynomials ipk{y) given in ( 12.191) 
are calculated by using (14. 8p . 

(4.9) A fc = -co(k + 2)k...[k + 2 - 2(m - 2)}(k + N + 2)(k + N)...[k + N - 2{m - 2)} 

for k = 2 (to — 3), 2{m — 2), ... . Using the eigenfunction expansion in terms of complete 
and closed subset of polynomials {ipk} partially justifies the asymptotic centre subspace 
analysis of the corresponding rescaled equation (I2.9p . which yields the same ODE (12.241) 
and hence the asymptotics (I2.26p . Here in the critical case (I1.25P we still have ^-j- = f3N 
with (3 given by (14. 3p . Finally, we arrive at the asymptotic pattern (ll.22p . where 4 is 
replaced by 2to. 

4.4. The general case n / 1. We do not have such a self-adjoint operator, but anyway, 
once F > in B\ is determined, we obtain the radial eigenfunction ipo for Ao = from 
the scaling symmetry group (12.281) (the exponent - is replaced by — ) of equation (I4.4p . 
We can also guarantee that (14.71) has compact resolvent provided that n > is not 
large, so A = is an isolated eigenvalue. The rest of the centre subspace behaviour via 
the expansion (12.331) remains unchanged and leads to similar logarithmically perturbed 
asymptotic patterns. A rigorous justification is a hard open problem. 
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5. Logarithmically perturbed patterns in the Cauchy problem 



The asymptotic behaviour and similarity solutions for the TFE (11.121) or (11.241) posed 
in the whole space M. N x M + are less studied in thee literature. For n G (0, |), in the 
Cauchy problem (CP), the solutions exhibiting the "maximal regularity" at the interfaces 
are oscillatory and of changing sign. See [TTl [18] and the book [24"l Ch. 1] for correct 
meaning of the CP for thin film equations and further examples. For such solutions, we 
need to assume that u n in (14.11) is replaced by \u\ n . Therefore from now on in all the 
expressions and equations we use the convention that 

u n , f n ,v n ,w n , ... are replaced by \u\ n , \f\ n , \v \ n , \w\ n , ... and, 
u p , f p ,v p ,w p , ... are replaced by \u\ p ~ 1 u,\f\ p ~ 1 f,\v\ p ~ 1 v,\w\ p ~ 1 w,.... 

We must admit that solutions of changing sign are less relevant for many known physical 
applications of TFEs. Nevertheless, for general PDE theory, it is key and of principal im- 
portance to include the Cauchy problem and to show that the basic techniques developed 
above apply to these much more complicated oscillatory solutions. 

The idea of sign changing solutions of TFEs is straightforward. Indeed, the oscillatory 
properties of such solutions are a manifestation of the fact that TFEs (14.21) are "homo- 
topic", i.e., can be continuously deformed (e.g., as n — > 0) via non-singular uniformly 
parabolic PDEs with analytic coefficients (see details in [THl Sect. 14]) to the linear poly- 
harmonic equation 

(5.2) u t = (-l) m+1 A m u in R N x R + . 

By classical parabolic theory (see e.g. Eidel'man [14J), given initial data uq G L 1 , there 
exists the unique solution of the Cauchy problem for (15.21) defined by the convolution 

N 1 

(5.3) u(x, t) — b(x, t) * w , b(x,t) — t~*^F(y), y = x/t^, 

where b(x, t) is the fundamental solution of the operator D t — (— l) m+1 A m . For any m > 2, 
the rescaled kernel F = F(\y\) is oscillatory as y —>■ oo, so this property of changing sign 
is inherited by L 1 solutions of (I5.2p . Assuming a continuous (homotopic) deformation of 
a class of solutions of (11.121) as n — > + , this confirms that the TFE admits oscillatory 
solutions of changing sign at least for not that large n > 0. Continuity and homotopy 
concepts are effective for treating the Cauchy problem for higher-order TFEs; see other 
examples in |18]. 

Then the source-type solutions of the TFE take the same form (I1.19p . where the radial 
function F of changing sign solves the ODE (ll.20p with the convention (15. ip . We begin 
with the linear case n = 0, which by continuity is going to describe some properties of 
source- type solutions for sufficiently small n > 0. 

5.1. Properties of the rescaled fundamental solution for n = 0. The linear ODE 



(5.4) 



A(F) = —A 2 F + \ VF ■ y + lj F = in K 
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is precisely the elliptic equation for the rescaled kernel F of the fundamental solution in 
(15. 3p . Therefore the similarity profile F(y) exists and is unique under the assumption 

(5.5) jF(y)dy=l 

(in view of existence- uniqueness of the fundamental solution). 

Let us next describe an important relation between similarity profiles for the FBP and 
the Cauchy problem. Without loss of generality, we consider the case N — 1, where on 
integration once (15 .4p takes the form 

(5.6) F"' = \Fy. 

It is easy to find all decaying profiles corresponding to the CP with the exponential WKB J 
asymptotics as y — > +oo, 

(5.7) F(y) ~ y~^e ay4/3 , with a satisfying a 3 = \ (|) 3 . 
There exist two complex conjugate roots for exponentially decaying profiles 

(5.8) a± = -§ 4-3(1 iiv 7 ^) = -c x ± ic 2 . 

This yields a two- dimensional bundle of oscillatory solutions with the behaviour 

(5.9) F(y) ~ y~^e~ Ciyi/3 [A 1 cos(c 2 y^) + A 2 sin(c 2 y^)] as y — > oo, 

where A\ and A 2 are arbitrary constants. The algebraic factor y^ 1 ^ 3 is obtained by a 
standard asymptotic WKBJ method. We observe here the periodic behaviour with a 
single fundamental frequency (a result we will refer to in the TFE analysis below). 

Proposition 5.1. For N = 1, the rescaled profile of the Cauchy problem F = F^ given 
by (5.4), (I5.5|) is the limit of FBP similarity profiles on bounded intervals, 

(5.10) F OQ =\imF k , 
where each F k (y) is defined on interval (—yk,Vk)> 

(5.11) F k (±y k ) = F' k {±y k ) = 0, and 

(5.12) y k = (iL A;) 1 (l + o(l)) as k - oo. 

Proof. The geometric aspect of such a property is obvious in view of the oscillatory 
behaviour in (15. 9p . The convergence as k — > oo follows from straightforward computations 
related to the whole exponential bundle including (15.9p and the growing counterpart 

F(y) = y~h aoy4/3 + ... , with a = §4~i. 

Then solving the FBP problem (15.111) yields the asymptotic equality cos (c2y k Z + const.) = 
0, whence the asymptotics (15.121) . □ 

We also expect the following Sturm property be valid: 

(5.13) F k (y) has precisely k zeros on (0,^). 

Such a zero- number property is easily seen for k ^> 1, but is not obvious for smaller k's. 
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FIGURE 2. The oscillatory CP profiles satisfying (15.141) . Parameters of 
shooting are F"(0) = -0.3379890 (n = 0), -0.3414702 (n = 0.2), 
-0.3490986 {n = 0.5), -0.3697143 (n = 1), and -0.4052680 (n = 1.5). 

5.2. Similarity profiles for n > 0: existence and uniqueness. 

Proposition 5.2. For N = 1 andn G (0, 1), i/ie (I5.il) m R admits a unique 

solution F G C 3 o/ umY mass. The solution F(y) is symmetric, compactly supported and 
is oscillatory near finite interfaces at y = ±y . 

Proof. For N = 1 the ODE (TOO]) has the form 

(5.14) \F\ n F"' = /3Fy, j/Gl. 
Dividing by |-F| n and setting l-Fl - "^ = g yields 

(5.15) (\g\ a g)"' = Pgy, v e R, « = ^. 

Then existence and uniqueness of a compactly supported solution F 6 C 3 for any n e 
(0, 1) follows from the results in Bernis-McLeod [6]. □ 

For n G [1, |) solutions of (I5.14p are less regular (see below), so the techniques in [6] do 
not apply directly, but we expect that the existence-uniqueness result remains valid and 
can be extended further to some interval n G see below. 

In Figure [2] we have shown these similarity profiles for some n > including the linear 
case n = leading to the ODE (15.61) for the fundamental rescaled profile. Here we observe 
convergence of the fundamental profiles as n — > + , which is justified rigorously if all the 
zeros are "transversal" and isolated except the last one; see below. 

5.3. Oscillatory properties via periodic orbits. We next describe the oscillatory 
properties of such changing sign profiles F(y) near interfaces. We rescale F to have that 

supp F = [—1, 1]. 

18 



It was shown in [T7] that the asymptotic behaviour of F(y) satisfying (I5.14p near the 
interface point y — ► 1~ is given by the expansion 

(5.16) %) = (1-I/Ms), s = ln(l-j/), /i = f, 

where, after scaling /3«0, the oscillatory component satisfies the following au- 
tonomous ODE (we omit exponentially small terms): 

(5.17) 0"' + 3(/x - 1)0" + (3/x 2 - 6/i + 2)0' + n(ji - l)(ji - 2)0 + ^ = 0. 

Oscillatory periodic orbits: existence. We are now interested in periodic solutions 
0*(s) of f)5.17p . which according to (15.161) . can determine the simplest typical (and possibly 
stable and generic) oscillatory behaviour of solutions near interfaces when s = ln(l — y) — > 
— oo as y — > 1~. There are several classic methods of ODE theory for establishing existence 
and multiplicity of periodic solutions of finite-dimensional dynamical systems. These are 
various topological techniques, such as rotations of vector fields, index, and degree theory; 
see [32j Sect. 13, 14]. Another approach is based on branching theory, [121 Ch. 6]. In 
our case, such an n-branching approach is especially effective since for n = the unique 
solution F is the rescaled kernel of the fundamental solution (a rigorous justification of 
some aspects of branching for such degenerate equations can be a hard problem). We also 
mention papers [HH EH EI] containing further related references and methods concerning 
modern theory of periodic solutions of higher-order nonlinear ODEs. In general, equations 
like (I5.17P are a difficult object to study, and especially the main difficulty is proving 
uniqueness of such periodic orbits. Therefore, later on, together with analytic techniques, 
we will need also to rely on careful numerical evidence on existence, uniqueness, and 
stability of periodic solutions. 

It is curious that for n = 1, the unique periodic solution can be detected by a direct 
algebraic approach; see [T71 Sect. 7.4]: 



Proposition 5.3. For n = 1, the ODE (5.17\) has a unique T '-periodic solution, with 



(5.18) T = -21ns > 9 = 1.9248... , 

where 9 = 0.381966... is the unique root on the interval (0, 1) of the cubic equation 

(5.19) 9 3 - 29 2 - 29 + 1 = 0. 

Indeed, for n = 1, the nonlinearity in (I5.17P is sign0 and the ODE is linear in the 
positivity and negativity domain of solutions, 

0'" + 60" + 110' + 60± 1 = 0, 

so can be solved explicitly. Matching positive and negative branches leads to the result. 
Let us now state the main result concerning periodic orbits of the ODE (15.171) . 



Theorem 5.4. The ODE {5.11) admits a nontrivial stable periodic solution 0*(s) of 
changing sign for all 

(5.20) < n < n h G (f,n+), where n + = = 1.9019238... . 
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m=2, n=3/2: to gobal stability of periodic orbit 
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FIGURE 3. Convergence to the stable periodic solutions of (|5.17p for n = | for 
various Cauchy data posed at s = 0. 

Uniqueness of such periodic 0*(s) in the interval ( I5.20p is still open. 
Proof. For the interval 

(5.21) 0<n<§, 

the proof of existence is performed in [T7J p. 292] by a shooting argument. Numerical 
representation of periodic solutions is given therein on p. 294; see also [24, p. 143]. We 
need to point out the main two ingredients of the proof in [T7] : 

(i) it is shown that for exponents (15.211) no orbits of the dynamical system (DS) (15.171) 
are attracted to infinity as s —* +oo, i.e., all orbits stay uniformly bounded; and 

(ii) as a consequence, then (I5.17P is a dissipative DS having a bounded absorbing set. 
Dissipative DSs are known to admit periodic solutions in rather general setting [321 

Sect. 39] provided these are non-autonomous (so the period is fixed). For the autonomous 
system ( 15.171) . the proof in [T71 Sect. 7.1] was completed by shooting. Note that, in view 
of the last term, (I5.17P is not a smooth dynamical system and solutions are not locally 
C 3 -smooth. Nevertheless, as shows local analysis [TTJ p. 291], at least for n 6 (0,2), 
the nonlinearity is integrable to guarantee local extensions of solutions through generic 
"transversal" zeros. This means that the equivalent integral equation is well-posed and is 
composed from compact operators in a certain topology (this is necessary for application 
of classic methods of branching in Banach spaces, [4"2"l Ch. 7]). We continue to deal with 
the differential equation, where the justification of calculus is done by local analysis. 

It turns out that both properties (i) and (ii) also remain valid for n = |, so that a 
periodic solution 0* also exists and is stable; see Figure [31 For the extension of 0* to 
n > |, we will use the following crucial stability result: 
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Proposition 5.5. If the periodic solution 0*(s) of (5.17) persists for all | < < 3, then 
it is stable and hyperbolic on this interval. 

Proof. Note that, for n G (§,3), there exist two unstable constant equilibria of (I5.17P 
(5-22) ^ = ±[- ^-4- 2 ) ]" for "6 (|,3), 

and we expect a stable periodic motion in between. Consider the eigenvalue problem for 
the ODE (I5.17P linearized about the T-periodic solution 0* by setting = 0* + Y, 

Y'" + 3(/i - l)Y" + (3^ 2 - 6fi + 2)r + - - 2)Y + (1 - n)|0,|- n F = AY. 

As usual, assuming that A G C, multiplying this by the complex conjugate Y in L 2 (0,T), 
taking the conjugate and multiplying by Y, and summing up both yields 

(5.23) -30* - 1) / \Yf + //(// - 1)0* - 2) / \Y\ 2 + (1 - n) J |0*MY| 2 = *±* / |F| 2 . 

Since all the three terms on the left-hand side of (15.231) are negative for any | < n < 3, 
the result follows. The case n = | is similar since just the second term vanishes. □ 

Thus, by classic branching theory, [121 Ch. 6], stable hyperbolic periodic solutions are 
locally extensible relative the parameter n > §. In particular, using the hyperbolicity of 
0* for n = §, we conclude that the periodic solution exists in an interval n G [|, | + 6) 
with some 5 > 0, and the interval of existence must be open from the right-hand side. 

Finally, let us justify the estimate in (15.201) . To this end, we multiply (15.171) by 0^ and 
integrate over (0, T) to get for any n G (0, 2) 

-J(0:) 2 + (3/i 2 -6// + 2)/(0;) 2 = O, 

so that one needs 

This completes the proof of Theorem 15.41 □ 

On heteroclinic bifurcation. Since the periodic orbit 0*(s) remains stable and hyper- 
bolic in the whole interval of existence (I5.20p . the end point n = cannot be any kind 
of subcritical saddle-node bifurcation, at which two branches meet each other. Classic 
bifurcation and branching theory [32, H2] then suggests that at n = % the DS (15.171) 
undergoes a heteroclinic bifurcation when the period increases without bound (this claim 
needs further study and a full analytical justification); see standard scenarios in Perko 
[lOj Ch. 4]. Note that, by Proposition I5.5[ the heteroclinic orbit occurred remains stable 
and hyperbolic. 

Numerically, is given by 

(5.24) n h = 1.7598665026... . 

Figure H] shows formation of the heteroclinic orbit in both limits: as n — > (a) and n — > 
(b). This bifurcation exponent plays the important role and shows the parameter 
range of n's, for which many ODE profiles near interfaces are oscillatory except those 
that approach the interface point s = — oo the stable manifold of the constant equilibrium 
(15.221) . In the interval (15.211) . this manifold of orbits of constant sign is empty, so that all 
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the orbits near s = — oo are oscillatory and coincide with the periodic one 0*(s + so), where 
so G M is a parameter of shifting. Indeed, this also characterizes important oscillatory 
features of the PDE. Note that some kind of a "heteroclinic bifurcation" phenomenon also 
exists for the sixth-order (m = 3) and higher-order TFEs with more difficult mathematics 
involved; see pH Sect. 13] and HQ p. 142-147]. 

On ID shooting for n G (l,nh). As a key application of the above oscillation analysis, 
we have that according to (15.16p . for all n G (0, n h ), there exists a ID bundle of oscillatory 
orbits of changing sign 

(5.25) F{y) = (1 - #0*(ln(l -y) + s ) + ... , 

where Sq G R is an arbitrary parameter of phase shift in the periodic orbit 0*(s). Recall 
that, for the ODE (15.141) . we need to shoot just a single symmetry condition at the origin, 

(5.26) F'(0)=0 (F(0)^0), 

so the ID bundle (15.251) is well-suited for this. In view of oscillatory character of the 
behaviour in (I5.25p . it is not a great deal to prove the existence of such a Sq to satisfy 
(15.261) . while uniqueness (as expected) remains open. 

Further comments about n^. For any n > n^, the behaviour in the ODE (15.171) 
becomes exponentially unstable and we did not observe oscillatory or changing sign pat- 
terns. This suggests that precisely above n = n^, the ODE (and the corresponding PDE) 
loses its natural similarities with the linear one for n = (though a continuous homo- 
topic connection is expected to be still available, i.e., some local properties of solutions 
dramatically change at n-^). 

Thus, in the range n G (§, 3), (15.171) possesses the positive constant solution </>(s) = + 
given in ( 15.221) . This gives the behaviour ( 12.351) . so that, for such solutions, formally, the 
FBP and the CP may coincide in the ODE setting. But this is not the case for all the 
solutions since for n G (§,?i-h) there are other oscillatory profiles with a similar (actually, 
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a bit less) regularity at the interfaces, so that the CP demands oscillatory solutions, while 
the FBP can admit positive solutions; see more details in [T?| Sect. 9]. In the parameter 
range n > nh, the oscillatory behaviour is no longer generic, so we expect a certain 
improvement of the positivity preserving properties of the TFE, where the CP and the 
FBP may coincide; see further discussion in [T71 Sect. 9.4]. 

5.4. The TFE with critical absorption p = p . The formal asymptotics for the TFE 
( 11.51) . ( 11.61) is now calculated similarly using the centre subspace spanned by the eigen- 
function (12.291) . Of course, we then do not gain any explicit mathematics or symmetric 
operators as for n — 1 in the case of the FBP. 

The main ideas of the analysis can be extended to the 2mth-order case, where many 
aspects of source-type and general solutions of the Cauchy problem for the TFEs remain 
mathematically open. The oscillatory character of solutions near the interface for m = 3 
was studied in [HI Sect. 13]; see also [2U Sect. 3.7] for further examples for m > 3 and 
other oscillatory PDEs. 

5.5. Supercritical range p > pq. We use the same scaling (13.11) and obtain the exponen- 
tially perturbed rescaled PDE (13. 2j) . which suggests that the solutions behave as t — > oo 
as the source-type solution with a finite positive mass attained at r = +oo (no proof is 
still available). 

Acknowledgements. The authors would like to thank J.D. Evans for discussions on 
thin film models with non-conservative aspects, and A. Leger for efficient consulting the 
authors with numerical methods for higher-order ODEs. 



Appendix A. The linearised operator is not symmetric when n/1 

We prove that, in the FBP setting, the linearised operator (|2.1ip admits a self-adjoint exten- 
sion only when n = 1. Without loss of generality we consider the one-dimensional case, and we 
formulate first the following results we are already familiar with. 

Proposition A.l. The linearised operator \2. 11\ in M is symmetric in some weighted space L? p 
when n = 1. 

Proof. For N = n = 1, the linearised operator is given by 



(A.l) 



A' 



(f)Y = -(/¥"')' -(Yf'"y + l(Yyy. 



For this to be symmetric in L 2 p with some weight p > 0, we require that |36|, Sect. 1] 



(A.2) 



A' 



V)Y = ± l( P0 Y")"-(PiYy+P2Y] 
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Expanding the right hand sides of these equations and comparing coefficients yields the following 



system: 








(A.3) 


Y»" ■ 


-f 


_ Po 

p ' 


(A.4) 


Y'" ■ 


-f 


_ M 
p ' 


(A.5) 


Y" : 


= 


Po'-Pi 

p ' 


(A.6) 


Y' : 


-r + b = - 


(A.7) 


Y : 


fllll , 1 P2 

J + 5" p' 



We know the exact solution of the ODE for / when n = 1 (see (|1.2ip ): 
(A.8) f(y) = (a 2 - y 2 ) 2 for y G (-a, a). 

Substituting this into equation (|A.7|) yields p2 = 0. Equation (|A.6P yields p\ = C where C is 
a constant. Equations (|A.3|) and ()A.4p yield Pq = f and p = — f^ 1 ^ 2 - Equation (|A.5P is thus 
the consistency condition and is satisfied since it yields p\ = C (since p'^ = p\ = C). Thus the 
linearised operator for the thin film equation is symmetric if n = 1. □ 

Theorem A. 2. For N = 1 and 1, operator (\2. 1 ip is not symmetric in L 2 for any weight 
p > 0. 

Proof. The ODE for / > for any n > is 

(A.9) _ (rry + _i_( /y y = 0. 

The linearised operator (|2,lip is given by 

(A.10) A'(/)y = -(/T'")' - n{r- x Yf»)' + ^ (Yy)'. 

For this to be symmetric, we require identity (|A.2p to hold. Comparing coefficients yields 



(A.ll) 


Y"" ■ 


fn po 

J p ; 




(A.12) 


Y'" ■ 






(A.13) 


Y" : 


. n - 
u p ' 




(A.14) 


Y' : 


■■ -nf n - l r" + ^u 




(A.15) 


Y : 


-n(n-l)r- 2 f" 


-nf n - 1 f"" + ^. 


From this 








pI = r, 


Pi = Po, P = 


---r n/ \ p2 = p[- 


n{n - l)f n - 2 f" - nf n 


and the consistency condition 


is 





i+4 J 



(A.16) /*(/*r = -nr-V w + A». 

To see if this coincides with equation (|A.9P for some / we use a Taylor expansion of f{y) and 
check if (fA~T6|) and fAT9|) prod uce the same coefficients for /. To do this we set /(0) = 1, /'(0) = 
f"'(0) = and /"(0) = b G M\{0}, differentiate equations (|A.16P and (|A.9P the required number 
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of times and set y = 0. The expansions coincide up to the coefficient of y 3 but the coefficients 
of y 4 only coincide if 



(A.17) ^ ^-"t^XT-^ 

Since we require b S R \ {0} we must have n E (—4, 0) U (|, 2). This gives us a range of values 
of n, for which the linearised operator may be symmetric. To check whether it is we examine 
the coefficient of y 6 for (|A.16|) and (|A.9p . If the operator is symmetric, then the same value of 
b should be obtained as in the coefficients of y 4 for both equations. For the coefficients of y 6 to 
coincide we require 

/ A 10 s , n , , 2 v / 2V rl ( 9n3 - 40n2 - 188n+464) (3n-2) 
( A - 18 ) b = °> OT b = ± V 9n<t-40n«-188n*+464n > 

and since we require b £ R \ {0} we discard 6 = 0. For this b to coincide with (|A,17p we 
require n = |. This contradicts the fact that we must have n G (—4, 0) U (|, 2) for the linearised 
operator (|2.1ip to have a chance of being symmetric and admit a suitable (Friedrichs) self-adjoint 
extension. Hence the linearised operator is not symmetric if n ^ 1. □ 
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